Bose-Hubbard model with attractive interactions 
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We consider tiie Bose-Hubbard model of atoms in an optical lattice potential when the atom- 
atom interactions are attractive. If the lowest energy lattice sites are degenerate (such as in the 
homogeneous case), then, at a critical value of the interaction strength, a phase-coherent conden- 
sate becomes unstable to a quantum superposition such that the number distribution of each of the 
degenerate sites becomes double peaked. In the limit when the interaction dominates, the superpo- 
sition becomes macroscopic and has the form \tp} oc e"'''St'^|vac), where A'^ is the total number 
of atoms and the sum ranges over the energy-degenerate sites. 
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An optical lattice potential loaded with a Bose- 
Einstein condensate of neutral atoms was predicted by 
Jaksch et al.^ and recently confirmed by Greiner et aZ.Q 
to be a realization of the Bose-Hubbard model of con- 
densed matter physics 0. In addition to exemplifying a 
quantum phase transition, this system is ideal for creat- 
ing and controlling the quantum states of the atoms at 
the sites of the lattice potential^ The Bose-Hubbard 
model is described by a Hamiltonian of the formal 

i?=-J^6]6, + ^e,n, + ic/n,(n,-l) (1) 

where hi and fii = blbi are the annihilation and number 
operators of the mode localized at the ith lattice site. J 
is the hopping matrix to neighboring sites, oc Os is the 
strength of the on-site interactions due to s-wave scatter- 
ing and Ci is an energy off-set due to an additional confin- 
ing magnetic trap. In the large hopping regime, J 3> \U\, 
the ground state of this system is well described by a 
phase-coherent condensate and the number fluctuations 
are approximately Poissonian: Srii = ^ (nf) — iif « \Aj7' 
where — {fii). For U > 0, increasing the inter- 
action energy reduces the on-site number fluctuations, 
{Srii < Y^) and, at a critical value of U/ J, the system 
undergoes a Mott-insulator phase transition {Sui — > 0) 
and all phase coherence vanishes 2] . In the experimental 
realization, the simple global phase-transition picture is 
complicated by the presence of the confining trap which 
produces local Mott domains^- 

The use of neutral atoms opens up the possibility of 
exploring the Bose-Hubbard model with attractive inter- 
actions {U < 0) as certain species of atom interact via 
a negative s-wave scattering length and many atoms can 
be made to interact via a negative scattering length by 
using the technique of Feshbach resonance to alter the 
interaction potential^. Attractive interactions are par- 
ticularly interesting because they lead to an instability of 
a phase-coherent condensate 0, 183 .9, 10, 11, 12J. The aim 
of this letter is to describe the nature of this instability 
for the Bose-Hubbard model and it's relationship to the 



formation of novel quantum states. 

Certain aspects of the behavior of the attractive Bose- 
Hubbard model can be deduced from the results for a 
simple two-mode model of a condensate in a double-well 
potentially 0|. In contrast to the repulsive case, it 
is predicted that the number fluctuations will increase 
{Srii > y/rTi) as the magnitude of the interaction energy 
is increased [T^. The behavior as the magnitude of the 
interaction energy is further increased depends crucially 
on the single-particle energies of the two wells, ei and 
£2. In the case where the wells are asymmetric, ei 7^ £2, 
the energy is minimized by all atoms accumulating in the 
lower energy site. In the case of symmetric wells, ei = £2, 
Cirac et al. 14j and Steel et al. |13| have shown that the 
system is unable to choose which site to accumulate in 
and will form a quantum superposition of the two possi- 
bilities. This superposition state is associated with very 
large number fluctuations: Srii Ui. In the multi-well 
case considered here, we expect this superposition state 
to form between all sites that are degenerate with the 
lowest energy site. 

These results show that if the site energies are degen- 
erate, then a phase-coherent condensate will become un- 
stable to a superposition state. On the other hand, in 
the absence of a lattice potential, a condensate with at- 
tractive interactions is known to become unstable as the 
interaction strength is increasedll i, H IH IH El . It 
is convenient to consider two distinct types of instabil- 
ity: (I) a global implosion of the condensate wavefunction 
confined in a harmonic trap[3, and (II) local instabil- 
ities of an unconfined condensate accompanied by large 
density fluctuations 9, T^l . In the present case we find no 
evidence of an instability of type I, even with the addi- 
tion of a harmonic confining trap, and conclude that the 
Bose-Hubbard model becomes invalid before this type of 
instability can occur. However, the instability where a 
condensate gives way to a superposition state described 
here shows many similarities to an instability of type II. 

We can analyze the stability of a condensate in a lattice 
potential via a Bogofiubov type treatment [isl IT^ : in the 
limit of large hopping, J \U\, the ground state is well 
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approximated by a phase-coherent condensate described 
by a mean-field, and we can consider small fluctuations 
about this mean-field by making the replacement bi{t) = 
^-ifit/h^p, -t- 5i{t)] in the Heisenberg equations of motion 
for bi and neglecting all terms except those linear in the 
fluctuations (5,;. This results in the equations 

= -J^/3j + (e,-M + t/|AnA (2) 
(i,i> 

^ = - Jj] <5, + {2Um' + u- fifS, + UpfSj (3) 

for the mean field and fluctuations, respectively. Af- 
ter solving Eq.Q for (3i and /i, Eq.lPJ can be 
solved by making the Bogoliubov transformation 6i — 
J2k "i,fc6~"^**4+Wj*fce"^'=*(5]!, with the normalization con- 
dition J2i — Iwi.fcP — 1, and solving the resulting 
equations for ujk, Ui^k and Vi^k- 

Assuming a homogeneous one-dimensional lattice of M 
sites (with periodic boundary conditions) containing N 
atoms, equations Q and ^ can be solved analytically 
such that fi^Un- 2J, where n = \f3i\^ = N/M 
is the mean-field solution and the quasi-particle energies 
have the form 

hcjk = V£fe(efe + 2J771) (4) 

where Sk = 4Jsin^(afc/2). Here k — {2n/aM)m for 
m = -M/2, • • • , M/2. This shows that as U be- 
comes increasingly negative, the first quasi-particle en- 
ergy (coresponding to fc = 27r/aM) drops to zero at 
Un = — 2 Jsin^(7r/Af) and then becomes imaginary, sig- 
naling a critical point beyond which the lattice system 
is unable to support a condensate. An analytical ex- 
pression for the on-site number fluctuations can also be 
determined from this treatment 

[13 (see also |i3) as: 
M ^ hcuk 

k 

It is evident that at the critical point the number fluctu- 
ations l(SJl diverge. Comparing with Ref. ^ we see that a 
type II instability is formally very similar to the present 
case. 

In order to treat the strong interaction regime beyond 
the instability (where the Bogoliubov treatment breaks 
down), we have numerically calculated the exact ground 
state. In the Fock state basis the state space of the sys- 
tem is large: {N + M-iy./[Nl{M-l)\]. ButtheHamil- 
tonian is a very sparse matrix and so for small num- 
bers of atoms and sites we can calculate the lowest few 
eigenvalues and eigenstates by the Lanczos methodls^. 
The results of these calculations are presented in FigQ 
After the critical point of the Bogoliubov treatment the 
exact calculations show that the number distribution be- 
comes double peaked (corresponding to the formation 
of a superposition state) which gives rise to the sud- 
den increase in the number fluctuations shown in Fig^ 
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FIG. 1: Ground state of a homogeneous one-dimensional 
lattice. This figure shows Srii (solid line), Ci (dashed line) 
and Cjv (dash-dotted line) as a function of U/J. The inset 
shows the energy difference between the first excited state and 
the ground state. The dotted fines correspond to a Bogoliubov 
treatment. 



The two peaks of the distribution move further apart 
and narrow as the interaction is increased which reduces 
the single-particle correlation between neighboring sites, 
Ci ~ (blbi+i), but increases the A^-particle correlation, 
Cn = M(5j^&^i)/iV!. Finally, in the strong attrac- 
tive interaction limit, the results conflrm that the ground 
state is a macroscopic superposition of the form: 

1 

IV') - e^*4t^|vac), (6) 

and the number fluctuations become Srii — 
N^/M - l/M or w N / \/M for many sites. 

In the current experiments^ Q, an inhomogeneity is 
introduced to the lattice system by an harmonic magnetic 
trap which is used to conflne the atoms in space. In 
the one-dimensional case considered here, this gives rise 
to the single-particle energies = \[i — \{M + A)]^, 
where i = 1 to M, A is a measure of the curvature and 
< A < 1 is the offset of the lattice from the center 
of the confining potential. Figure |2] shows the results 
of exact calculations in the case of non-degenerate sites. 
Note that the number fiuctuations become large close 
to where the critical point would have occurred if the 
sites were degenerate (the system gets close to forming 
a superposition at this point) and then decrease as the 
atoms accumulate into a single site. 

We can also calculate the critical behavior of the inho- 
mogeneous system via a Bogoliubov treatment, which re- 
quires the numerical solution of the mean-field equations 
(O before solving the linear equations Q for the quasi- 
particles. In the non-degenerate case when < A < 1, 
for M odd, there is always one site, i = \{M + 1), with 
the lowest single-particle energy. In this case, no crit- 
ical point is seen and as the interaction becomes more 
attractive, the condensate-which has an approximately 
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FIG. 2: Ground state in the presence of a harmonic confining 
trap with A = O.IJ and A = 0.9. Figures (a)-(c) show rii 
(white bars) and Srii (black bars) at each site for various val- 
ues of U/J. Figure (d) shows 714 (dashed line) and Sn4 (solid 
line) as a function of U/J. The inset in (d) shows the en- 
ergy difference between the first excited state and the ground 
state. 
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FIG. 3: Minimum energy gap between the ground state and 
first excited state as a function of A for AI = 50 sites. The 
values of the curvature are A — 10^^ (solid line), A = 10~^ 
(dashed line) and A — 10~^ (dotted line). 



Gaussian spatial profile of width w-simply decreases in 
width until all atoms accumulate in site i — ^(M-f 1). In 
figure 13 we have plotted the minimum energy gap to the 
first excited state, Aii^niin, (found by varying U/J over 
a broad range of values) as a function of A. The finite 
value of this gap for A 7^ demonstrates the stability 
of the condensate in the non-degenerate case. We find 
no evidence of an instability of type I which would be 
expected to occur even in the non-degenerate case. 

Figure 0] shows the dependence of the critical point 
(where the first excitation energy becomes imaginary) 
on the curvature in the degenerate case: A = 0. The 
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FIG. 4: Plot of the critical value of U/J as a, function of 
the curvature of the confining trap. A, in the case of two 
degenerate sites at the center of the trap. The inset shows 
the width w of the condensate at the critical point. 



inset shows that for a broad range of trap curvatures, 
the critical point occurs when the condensate width is of 
the order of the site spacing. 

In a realistic situation it may be difficult (if not impos- 
sible) to create a lattice with exactly degenerate sites. 
However, in current experiments, the high interaction 
regime is reached by a dynamic process whereby the rel- 
evant parameter, K{t) = \U{t)\/J{t), is increased at a 
certain rate 7 (this is normally in the adiabatic regime 
so the system remains in the ground state) . In this case, 
close to K{t) ~ 1/A'' (where the critical point would be 
if the sites were degenerate), if one increases K{t) at a 
rate faster than the oscillation frequency between two 
different sites, 7 » (ei — £2)/?*, (but still slower than the 
tunneling rate) then, for short times, the system will be 
unaware of the inhomogeneity and it is possible to effec- 
tively "capture" the superposition before the atoms are 
able to tunnel into the lowest energy site. To illustrate 
this point we have numerically solved the Schrodinger 
equation for a two-mode model with the time-dependent 
interaction strength U{t) = Uf{l — e~^^). The results 
of this simulation are shown in Fig[Sl The large number 
fluctuations for a sufficiently large 7 (dotted line) con- 
firm the formation of a superposition even though the 
single-particle energies in the two wells differ. 

Experimental realization of the interaction dominated 
regime of the Bose-Hubbard model with a large number 
of atoms is complicated by the fact that the localization 
of all N atoms at a single lattice site may render the Bose- 
Hubbard model invalid unless the magnitude of the scat- 
tering length is small. Approximating the potential at 
each lattice site by a harmonic potential of length Oho, the 
interaction strength must satisfy niOho 3> {hi(ni — l))\as\^ 
such that the interaction does not alter the shape of the 
localized mode- functions at each site [In fact, we can con- 
tinue to use the Bose-Hubbard model beyond this in- 
equality with renormalized parameters 19j, the ultimate 
limit being the stability of the localized mode-function 
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FIG. 5: Non-adiabatic evolution of a two-mode model from 
the ground state a,t U(t — 0) = with the parameters 
Uf = -0.05J, £2 - ei = 0.0005J. This figure shows the 
time evolution of the number fluctuations, Sni for the rates 
7 < 0.005 J/?i (solid line) and 7 = O.OIJ/?i (dashed line). 
The inset shows the two-peaked structure (coresponding to 
a superposition state) of the number distribution, P{ni) — 
\{'ip{t)\ni, N — ni)f , at the point corresponding to the arrow. 



against a type I instability: Oho > JT^ilcisI]- A possible 
method to overcome this is to load a very deep optical 
lattice (so that J is small) with repulsively interacting 
atoms and use Feshbach resonance to slowly tune the 
interaction through zero so it becomes just slightly neg- 
ative, as described in Ref.[20l| for ^Li. This method has 
the additional advantage that it will minimize three-body 



loss of atoms which scales as alpH. 

In conclusion, the ground state of the attractive Bose- 
Hubbard model displays behavior fundamentally differ- 
ent from the repulsive case. In particular, if the low- 
est energy sites are degenerate then, at a critical value 
of the interaction strength, a phase-coherent condensate 
becomes unstable to a quantum superposition such that 
the number distribution at each degenerate site becomes 
double peaked. The atoms have a tendancy to accumu- 
late at a single site in order to minimize the interaction 
energy but they are unable to choose which site due to the 
energy degeneracy and so form a superposition of all the 
possibilities. Interestingly, our results suggest that the 
Bose-Hubbard model becomes invalid before an instabil- 
ity of type I (for a confined condensate) occurs and that 
an instability of type II (for an unconfined condensate) 
corresponds to the formation of superposition states due 
to the homogeneity of free space. In an experimental re- 
alization, atom loss (or absorption imaging) will destroy 
the superposition by tending to localize the atoms at one 
site. Superpositions such as © can be destroyed by the 
loss of just one atom which will "collapse" the quantum 
state of the atoms to one of the degenerate lattice sites. 
Less macroscopic superpositions, such as those formed 
just after the critical point, will be more robust against 
loss 14]. Methods of non-destructive detection of these 
superposition states will be the topic of future work. 

The authors would like to thank A. Kawaguchi, Y. 
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